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1. INTRODUCTION 

The IAU/IUGG Working Group on "Nutation for a non-rigid Earth", led by Veronique 
Dehant, won the European Descartes Prize in 2003, for its work developing a new model for 
the precession and the nutations of the Earth. This model (MHB2000, Mathews et al. 2002) 
was adopted by the IAU (International Astronomical Union) during the General Assembly in 
Manchester, in 2000. It is based (i) on some improvements for the precession model (with respect 
to the previous one of Lieske et al. 1977) owing to the VLBI technique (Very Long Baseline 
Interferometry) , and (ii) on a very accurate nutation model, close to the observations. With this 
prize, the Descartes nutation project could offer for international scientists some grants, to be 
used for further improvements of the precession-nutation Earth model. At the IGG (Institute of 
Geodesy and Geophysics), with the OCCAM 6.1 VLBI analysis software and the best data and 
models available, we re-analyzed the whole VLBI sessions available (from 1985 till 2005) solving 
for the Earth Orientation Parameters (EOP). In this paper we present the results obtained for 
the EOP and more particularly for the nutation series. We compare them with the other IVS 
(International VLBI Service) analysis centers results, as well as with the IVS combined EOP 
series from the analysis coordinator. The series are in good agreement, except for the polar 
motion coordinates that show a shift with respect to the other ones and that we discuss here. 
Finally, we analyse the nutation series in the framework of the free core nutation (FCN) effect. 



2. COMPUTATION 

In this study we use 2944 VLBI sessions from 1982 till 2004, described as suitable for de- 
termining the EOP by the IVS Analysis Coordinator. They are computed at the IGG, from 
the NGS-format files corrected on the basis of the ECMWF meteorological data, in order to 
improve the height component of the stations and the determination of the tropospheric pa- 
rameters. Clock breaks and reference clocks were investigated and taken into account for each 
session (Heinkelmann et al., 2005). The VLBI analysis software used at the IGG is OCCAM 
6.1, in which the classical least-squares method based on the Gauss-Markov model (Koch, 1997) 
is implemented. We used (i) the terrestrial reference frame ITRF2000, (ii) the Vienna Mapping 
Function (VMF; Boehm & Schuh 2004), and (iii) a cut-off elevation angle for the troposphere 
set to 5 degrees. We estimated atmospheric gradients, clock parameters and zenith delays, as 
the same time we solve for the 5 EOP: one per session, and its rate for the pole and UT1. But 
we do not estimate the stations and sources position. 



3. DISCUSSION 

In Figure [TJ we show the results obtained for the polar motion (x p component) compared to 
the IVS combined rapid solution (series ivs03rle). This series was obtained using the IVS-R1 
and R4 sessions, very good ones for determining the EOP, occuring after the year 2000. We 
can notice a shift in our series with respect to the IVS combined one. We wondered if this 
shift could be solved removing Tigoconception station. And we can realise that it is in better 
agreement after removing it (see Fig. []]). This can be explained by the fact that we used 
ITRF2000 terrestrial reference frame, in which Tigoconception position is not well determined. 
But it still remains a difference at the end of the data time span. In Figure El we plotted 
the Fourier spectrum of the Celestial Pole Offsets (d^, de) we obtained, with respect to the 
IAU2000 nutation model. Depending on the data time span considered, we do not obtain the 
same periodical signals. The main ones are summarized in Table [TJ investigating the longest data 
time span possible (between 1982 and 2004). For further studies about the Free Core Nutation 
effect (FCN) , we will investigate a wavelet analysis on the prograde and retrograde parts of the 
Earth nutations. In the future, we will be able to use also the intensive VLBI sessions in our 
computations. 




Figure 1: Polar motion x p component, ob- 
tained with OCCAM 6.1 and compared to the 
IVS combined series. 



Figure 2: Fourier spectrum of the celestial pole 
offsets obtained with OCCAM 6.1, with re- 
spect to IAU2000 model. 
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Table 1: Amplitudes of the biggest periodical signals for (d^,de), with respect to IAU2000. 
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